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The main aim of this study is to support design and development processes of 
advanced fuzzy-logic-based controller for smart buildings e.g., heating, 
ventilation and air conditioning, heating, ventilation and air conditioning 
(HVAC) and indoor lighting control systems. Moreover, the proposed 
methodology can be used to assess systems energy and environmental 


performances, also compare energy usages of fuzzy control systems with the 
performances of conventional on/off and proportional integral derivative 
Keywords: controller (PID). The main objective and purpose of using fuzzy-logic-based 
model and control is to precisely control indoor thermal comfort e.g., 
$ temperature, humidity, air quality, air velocity, thermal comfort, and energy 
Fuzzy logic balance. Moreover, this article present and highlight mathematical models of 
Indoor thermal comfort indoor temperature and humidity transfer matrix, uncertainties of users’ 
PID controller comfort preference set-points and a fuzzy algorithm. 
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1. INTRODUCTION 

The integration of real-time events with indoor air temperature, humidity, the level of carbon dioxide 
(CO2) and indoor illuminance level are not a straightforward task using proportional integral derivative (PID) 
controllers. Thus, a multi-agent control system, are designed and reconfigured in this paper rather than using 
one multi controllers of each separate control signal. In addition, it is quite difficult to apply and design model 
free controller to real-time indoor environment control system based on surrounding real life event and real 
time thus modifying adaptability is included in this study in order to increase controller feasibility. 
Consequently, suitable mathematical models of controlled goals are essential throughout the design and testing 
phase, the proposed control strategies. Accordingly, this chapter discusses the mathematical models of the 
building’s indoor atmosphere including the model of indoor temperature, humidity, the concentration of CO2 
level and the indoor illuminance level. The chapter also discusses the mathematical models of lighting and 
heating, ventilation and air conditioning (HVAC) system and systems’ energy usage and savings. The recent 
control methods that extensively employed in indoor thermal comfort control mathematical models. The 
current research, a fuzzy-logic-based controller is designed and modelled in order to control indoor 
temperature, indoor humidity, indoor air quality and indoor illuminances by integrating real life events such as 
the usage of ambient air and natural light (day-lighting) into the existing artificial indoor environment while 
using minimum energy. 
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The mathematical models have been investigated heavily throughout literature such as [1] have 
modelled and experimentally investigate indoor thermal comfort. Also, same authors have highlighted indoor 
comfort parameters and its effect on indoor air quality [2]. Lately [2]—[6] have introduced necessary parameters 
in considering indoor thermal comfort in fuzzy logic controllers. However, the most of those research activities 
lack of accurate mathematical models. Accordingly, this article presents and highlight mathematical models of 
indoor temperature and humidity transfer matrix, uncertainties of users’ comfort preference set-points and a 
fuzzy algorithm and the integration of daylight with existing artificial light. 


2. INDOOR ENVIRONMENT MATHEMATICAL MODEL 

Here, the controller’s modelling representing indoor climate level containing building’s inner space 
temperature, building’s indoor humidity, building’s indoor air quality (LAQ) and comfort index is presented. 
The proposed control system which is based on add on multi-agent fuzzy based control system into existing 
PID controller algorithms are detailed and presented in the following sub-sections. 


2.1. Modelling of space indoor temperature 

Inner spaces’ temperature of a building is a very important factor in indoor environment monitoring. 
The models of in space air temperature is able to predict the daily change of indoor temperature with time. 
According to [7], the in-space temperature is influenced by the atmosphere temperatures in inside climatic 
conditions. The modelling process of inner air temperature is built using the regression analysis method by 
evaluating conduction, radiant and convection heat transfer. Accordingly, indoor temperature can be evaluated 
based on (1) taking into account energy conservation principle: 


dt 
Pa- Cp. Vi -- = Qhi a Uwâw(T; ga To) (1) 


where pa is air density in (kg / m°), Cp represents heat capacity of air in (J/kg-C°), Vi is room volume in (m°), Ti is 
inner doors temperature in C°, t is time in (s), Qn is HVAC work rate, To is outdoor temperature in C°, Aw (m7) 
is the area of walls, Uw is overall heat transfer coefficient in (W/m? C°). Uw can be calculated using (2). 


Uy = aa (2) 


ha kb hB 


Here K, represents the walls heat conductivity, units is (W/m - k), dw is walls thickness in (m), ha and hg are 
convection heat transfer coefficients on the side of each wall in (W/m? : k). Hence Laplace transform is used to 
transform (3) into (3) and then into (6) taking into account the convection heat transfer coefficient is constant [8]. 


1 


T(s).S = z h(S) (3) 
T(s).S = AT Qh1(S) + aa T(S) (4) 
UwAw 

meee OSIEN (5) 


Then (4) can be written as (6): 


Ti(S) E Kit 
QRh1(S)  Tit(S)+1 


(6) 
Here T; is the time constant and heat gain is a DC gain of the transfer function is K4. 


2.2. Modelling of inner space humidity 

The mathematical model of indoor humidity can be obtained in a similar way to the of indoor 
temperature modelling. The humidity of a conditioned space room is considered to be equal to the mean 
humidity that are provided by sensors. The humidity of a space depends on the humidity of the provided air 
from the inlet air and HVAC ducts, room temperature, and occupant’s number. Using energy conservation 
principles, indoor air humidity can be expressed using (7) [9]: 
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dhi 


Pa- ViH.— = Qh2 — fH palhi — ho) (7) 


where vapor enthalpy of vapor in H ( kJ/kg), fis the flow rate volume in (m3/s), hi is the inner space humidity 
in (g /kg), ho is the outdoor humidity in (g/kg) and Qh2; work rates in (W) and drt is time constant. 
The (8) is resulted using Laplace transform in order to solve the previous differential (7). 


1 
PaViH 


H,(S).S = 


Qh2(S) — Ž (H,(S) ~ Ho(S)) (8) 


Here, (9) is resulted from rearranging (8). 


1 
PaViH 


(s+ 4) H,(s) = Qh2(s) + 4(H,(S) a 


Taking into account that Lh, (S) = 0 as the change in indoor enthalpy is neglected due to its’ little effect 


and, then the transfer function is given in (10): 


H;(S) = K 
QRh2(S) Tin-(S+1) 


(10) 


1 
PaViH 


Vi. : . f P 
where Tin- T is the time-constant and kiņ= is the heat gain transfers function. 


2.3. Modelling of indoor thermal comfort conditions 

The modelling of indoor thermal comfort conditions and environment is defined by ASHRAE 55 and 
ISO 7730. In addition, comfort in an indoor environment, indices also the way of measurement methods are 
presented here. This section also shows most of parameters that affects local thermal comfort level such as air 
speed, thermal radiant gain, and vertical temperature difference as shown in Figure 1. There are a few research 
activities available concerning the mathematical model of thermal comfort [10], [11]. These authors have 
investigated the implementation of mathematical models in the control process of smart buildings and thus 
improving the thermal prediction level of pivot mean value (PMV) model in buildings. 


Temperature difference 


Local Thermal model 
Vertical 


Figure 1. Indoor standards model 


2.4. Mathematical model of thermal comfort 

Air speed is linked to air sensible heat gain caused by convection heat transfer and latent heat gain 
caused by evaporation and, thus, the sensation of thermal comfort is affected by air current. For example, 
occupants expose a more indoor air temperature the moment indoor air velocity is decreased. Accordingly, 
throughout hot season, ventilation is increased to cool the conditioned space rather than running the cooling 
system which means the cooling energy consumption is decreased. Throughout cold season, a higher air 
velocity is required for operating heating system which means heating energy consumption is higher. In order 
to control air velocity influence, air velocity limitation is applied as indoor air velocity ranging between 
0.6-0.90 m/s in hot season and under 0.150 m/s throughout the cold season taking into account temperature 
fluctuation. Simultaneously, in two different situations, occupants encounter the same heat loss and thus the 
feeling of higher temperature unsatisfactory where air temperature encounters high varying in respect to time, 
which is known as air turbulence [10], [11]. The pressure in duct percentage PD% can be obtained based 
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on (11) which defines the draft risk (DR) which expresses the percent of people dissatisfied due to draught. 
The implied set point temperature is in the range of 22—26 °C, Air speed velocity varies between 0.050 m/s 
and 0.40 m/s, and winds concentration is ranged from 0% to 70% [12]. 


DR = (34 — t) (v — 0.05)%62(.37v. Ty + 3.14) (11) 


3. THERMAL RADIATION AND VERTICAL TEMPERATURE DIFFERENCES MODEL 

Thermal radiation is defined as the heat energy gained by the mean of emission of electromagnetic 
waves. Asymmetric thermal radiation occurs when residents of an inner space environment are indorsed in to 
a heat source on one part of their body for substantial time, leading to a certain level of un-satisfactions or 
discomfort. The effect of current thermal radiation is obtained by two techniques. 

Measuring in 2 opposite directions and utilizing transducer to record radiation that affects a tinny 
point from the matching hemisphere and measuring surroundings surface temperature in indoor environments 
and evaluate the increase of radiant temperature. Consequently, (12) shows the percent of PID caused by warm 
ceiling [13]. 


100 


PD = 1pel28470.174Atpr) 


5.5 (12) 


In addition, the percentage of PD caused by a cold ceiling is given by (13). 


100 


1+e(993-0.54tpr) 


PD = (13) 


In a similar way, the percentage of PD caused by cold and hot walls are given in (14) and (15) respectively 
[13]: 


100 
ET . 
100 
BD epea p 


where Atpr is the increment of radiant temperature in °C. The vertical temperature differences affect PD for 
various increments of temperature, as shown in (16): 


100 


PD = 1+e(576-0.8564Atpr) 


(16) 


were At is the vertical temperature difference in °C. 


3.1. Thermal comfort index 

American Society of Heating, Refrigeration and Air Conditioning Standards (55) have outlined the 
influences of tinner space temperate and inner environment on residents’ satisfaction under hot and humid 
atmosphere also under misty atmosphere based to the model of humans’ metabolism enhanced by that human’s 
thermal behavior which is influenced by different environmental conditions and several sole conditions as there 
is a definite difference between the PMV model and the thermal sensation, as presented by (metabolic rate and 
clothing) [14]. The difference, compared to the neutral temperature (comfort level) is nearly 1.4 °C due to heat 
sensing is acquired from occupants’ survey located in certain environment taking into account that, PMV 
method utilize heat transfer as demonstrate in (17): 


PMV = at + bp — c (17) 


where a, b and c are coefficients representing occupants’ gender and their function of spent time as shown in 
Table 1. Meanwhile the error in evaluating the body neutral temperature is related to a PMV issue in 
determining the metabolic rate and the clothing factor (clo) which represents the insulation of the body [15]. 

Environments’ researches of Kansas’s State University helped developing a model of PMV that is 
able to evaluate inner thermal comforting an institutional building in order to quantify the parameters of PMV 
[15]. This research project developed different regressions and consequently, the thermal sensational value 
Tsens Which is presented in Table 2 as in (18): 
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Tsens = 0.305T + 0.99clo — 8.08 (18) 


where Tsens the sensational temperature in °C, T represents inner temperatures in °C also clo is clothing factor 
in clo. 


Table 1. Coefficients (a, b and c and gender of the resident) 


Time s / Genders a b Cc 

1 h/Males 0.2190 0.2330 5.670 
Females 0.2719 0.24890 7.240 

Both (male sand females) 0.2445 0.2490 6.470 
2 h/Males 0.2220 0.2690 6.020 
Females 0.2820 0.2990 7.690 

Both (male and females) 0.2510 0.2390 6.850 
3 h/Males 0.2130 0.2990 5.940 
Females 0.27560 0.2560 8.6219 


Table 2. Coefficient a-b-c and genders 


Thermal sensation Thermal sensations 


3.0 Warmth 
2.0 Hot 

1.0 Near hot 
0.0 Neutralth 
-1.0 Fresh 
-2.0 A bit mist cold 
-3.0 Coldth 


It is worth mentioning that the previous models that were endorsed to assess overall thermal 
comforting. Researches’ activities have shown that adaptive-model could be utilized to calculate sat thermal 
conditions inside the space atmosphere based on functions of outer ambient climates. This model is able to 
utilize where occupant are in resting mode or seated activity conditions also the Met level is ranging between 
1-1.3 Met. The ASHRAE proposed model is based on windows that are the medium of controlling thermal 
conditions as shown in (19): 


T, = 17.8 + 0.31T (19) 


here Te is comfort temperature in °C, T, is outer air temperature in °C and 7; is the mean inner air temperature 
in °C. Percentage of dissatisfied persons with local thermal comfort (PDP) is given by (20): 


_ 100 
~~ 44 e(-3.58+0.18(30-t)+0.14(42.5—0.01P)) 


PD (20) 


here P is number of persons. 


4. INNER TEMPERATURE’S AND HUMIDITY’S TRANSFER MATRIX 

Generally indoor environment temperature variation and indoor humidity variation are connected, it is 
important to evaluate the two variables jointly and thus establishing indoor temperature and humidity transfer 
matrix. In this research, the system output control variables are inner temperatures humidities. The output 
matrix is given by (21): 


Y = [T; - HJ" (21) 


where T; is the inner temperatures in °C and H; is the inner humidities in g/ kg. Considering the cooling system 
including humidifier supply power are the inputs while the output are the controlled temperatures and 
humidities. Accordingly, the input matrix given by (22): 


Y = [Qh; — Qh)” (22) 


here Qn: is the cooling system power supply in kWh and Qn2is the dehumidifier power supply in kWh. The 
model of the two input parameters and the two output parameters is shown in Figure 2 where T-controller is 
the temperature controller and H- controller is humidity controller. 
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T-controller |Q# B11 > 
j B21 


H-controller 37 
_>| > B22 
Rh p On Sa 


Figure 2. Inner temperatures and humidity’s single input, and output model 


R, Š 


In the previous model, the model transfer-matrix which evaluates the relationship between output Y 
and input X is given in (23): 


By a | Fa 
23 
a] Femme lie Sa 
where 7;, is the indoor temperature controller output in °C, Hic is the indoor humidity controller output in g/kg, 
Qnı is work rate of heater in W, Qn2 is work rate of humidifier in W and B’s are the matrix transfer function. 
Hence, the air supply flow rate is constant; the humidity is not affected by air cooling process power as it is 


influenced by the dehumidifier’s process power. In short, the mathematical description of the previous 
paragraph means that G22=0. Simplifying (26) will give (24). 
B11 By, Be 

HE [o B22 | Fa (24) 
4.1. Mathematical model of thermal comfort 

In this thesis, the indoor temperature and humidity is controlled by integrating outdoor ambient air 
into indoor once its temperature and humidity are suitable to use. In other words, if the outside air temperature 
is ranging from 20-25 °C and its humidity is ranging between 40%-70%, the ambient outdoor air can be 
pumped into the building, and, hence, less power is needed to power HVAC system. The indoor temperature 
and humidity are considered to be difficult as both temperature and humidity are linked together and 
consequently any indoor humidity changes, will endorse changes in the indoor temperature. Therefore, indoor 
temperature and humidity single loop control system is not able to cover the optimum control necessity and 
also it will produce poor control functioning. Accordingly, the design and implementation of decoupling 
control strategy will remove the link between indoor temperature and indoor humidity. Later, once the indoor 
temperature and humidity are not linked to each other, using single-loop controllers to control indoor 
temperature and humidity control can be satisfactory. Figure 3 represent a two input and two output control 
system where the input and output variables affect each other’s, in the figure R, and Rp are the system inputs 
(the outdoor air temperature and humidity). X; and Xn are the controllers’ temperature control outputs and 


humidity control outputs respectively while Y, and Yp are the system temperature and humidity output, 
respectively. The input matrix is obtained from (25) as shown in (25): 


Xtn a [Qn = Qnr2]” (25) 


where Qn; is the cooling system power supply in kWh and Qn2is the dehumidifier power supply in kWh. 
Furthermore, the output matrix is given by (26): 


Xen = [T; — H" (26) 
where T; is the indoor temperature in °C and H; is the indoor humidity in g/kg and hence the matrix transfer 


function can be given by (27) and (28): Adding the decoupling elements; F7;, E12, E2; and E22 as shown in 
Figure 4 to (30) will result a new transfer matrix as shown in (28): 
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“J [2 By | Ee] 
= 27 
Yp] lea, Balle On 
| = a fal Eo es K (28) 
Yn B21 B22 E21 E22 Xn 


Rt y- Xx: 
nhs, T-controller ix 
5 


Figure 4. Two-input two-output control system with decoupling 


Assuming that the decoupling elements; F7;, E12, E2; and E22 are equal to the inverse of matrix transfer 
function as in (29) and (30): 


Ey, 7 | ka A 

= 29 
En Eglo [Ba Bz 29) 
Y, _ By, O07 /X: 
l = fF Boa | G8) 


Consequently, the two-input and two-output control system is represented by two single-input single-output 
control system (SISO) as demonstrated in Figure 5. 

As mentioned previously, temperature and humidity are given by (24). and can be noted that B2);=0, 
so the flow chart of the control system can be summarized in Figure 6 which contains two control loop: 
temperature control loop and humidity control loop to ensure outdoor conditions are used by removing the 
effect of Bz;. Consequently, the decoupling control element of the system E72 is configured as in (31) which 
can be added to the previous control system as shown in Figure 7. 


E12(S) = O (31) 


B11(s) 
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It is may possible controlling inner temperature and inner humidifies utilizing 2 single loops controller that is 
capable to add outer ambient atmosphere in treated inner environment [16]—[19]. 


T-controller 


—+( H-controller 

Rh 
Figure 5. SISO Equivalent control system 
R On T; 
: T-controller +) Bl 
B12 + 

i—s( )——_s}H-controller B22 — Or 

Rh Qnm Hi 


Figure 6. Temperature and humidity control system 


R, 
T-controller 


— Or 


Figure 7. Temperature and humidity control system with decoupling strategy 


5. RESULT AND DISCUSSION 

System’s simulations medium evaluated HVAC’s System also the illumination systems energy usage 
knowing that residents’ occupancy changes. ADDON fuzzy PIDs’ controlling systems which readjusted so it 
could carry on all calculations. The fuzzies controlling input is (eR) also errors variation (4 eR) that highlights 
variations of buildings former residents and building’s todays residents. Additionally, memberships formations 
of ADDON (fuzzy PID) controlling inputs as well as outputs all show in Figure 8. The controlling input-output 
uses the stated functions value: negative-high (NH), negative-medium (NM), negative low (NL), zero (Z), 
positive low (PS), positive medium (PM) and positive high (PH) are demonstrated in Table 3. The output of 
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ADDON (fuzzy PID) controlling systems are basic energy variations depends on residents occupancy density 
that endorses keeping inner thermal-visual comforting level. The research, BMS-fuzzy logic controlling 
systems of institutional spaces under Australian Rockhampton campus of Central Queensland University have 
been evaluated, designed, and built. The performance of the proposed control system was numerically 
simulated and analyzed using MATLAB, DAYSIM and EnergyPlus software [20]—[22]. The study presented 
a comprehensive understanding and assessment of control strategies of smart buildings application using 
advanced fuzzy logic based control system. 

The study investigated existing conventional control systems including control for indoor 
temperature, humidity, air quality, HVAC system and the lighting system. Help designers to designed and 
developed models for fuzzy-logic based controllers of the reference building. The mathematical model of 
proposed fuzzy-logic-based control system was developed based on real-life events such as the utilizing 
un-conditioned ambient fresh air, day-lighting, and the ability of performing head count for indoor environment 
such as indoor temperature, humidity and thermal comfort [23]-[26]. In addition, the study will help in 
developing methodologies and strategies to achieve potential energy savings in the reference building taking 
into account of indoor environmental conditions. The existing and proposed control systems including 
conventional PI controller, fuzzy logic controllers and add on multi-fuzzy-logic-based controllers were 
simulated using energy plus, DAYSIM and MATLAB [27]. 


2 
E 
|= 
D 
§ 
= 
5 
Y 
D 
© 
& 


0 
eR/Delta eR 


Figure 8. Membership functions of various occupancy levels 


Table 3. Fuzzy control rules for local various occupancy levels 
Energy changes eR 

NE’s NM’s NLĽ’s Z’s PL’s PM’s _ PL’s 

(AeR) NH’s NH’s NR’s PHs PH’s PHs PHs PH’s 

NM’s NE’s NM’s Z’s PM’s PM’s PH’s_ PH’s 

NL’s  NH’s NM’s NE’s PHs PM’s PH’s PH’s 

Z NE’s NM’s NH’s Z’s PL’s PM’s PH’s 

PL NH’s NH’s NM’s NE’s PLĽ’s PM’s PH’s 

PM NE’s NEs NM’s NM’s ZEs PM’s PH’s 

PL NE’s NH’s NE’s NHE’s NH’s PHs PH’s 


6. CONCLUSION 

Occupants spent the most of their time indoor, thus indoor environment has major effect on occupants’ 
health and productivity. The main factors of indoor climate aspects are indoor temperature, humidity, inner 
illumines, and inner air quality. Nevertheless, almost inner climatic features are calibrated by utilizing regular 
controlling systems that its functions depends on ON/OFF controlling system and PID controlling system. As 
mentioned earlier, the conventional controlling system are not a suitable solution in saving energy as the 
functions of HV ACs’ systems and illuminations’ systems are non-linear. Therefore, the utilization of fuzzy 
controlling systems in smart buildings will be energy efficient that is able in saving energy and also 
accordingly, able greenhouse gas emission reduction and reducing its negative impact on the environment. The 
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structure of the newly proposed fuzzy logic based controllers contains two main segments; a market available 
conventional PID control which normally control indoor air temperature and an add-on fuzzy logic controller 
which is able to integrate the newly controlled real life events such as outdoor air temperature, humidity, and 
daylight usage to the PID controller based on current indoor situations. 
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